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ABSTRACT 

Context. The dust properties in dwarf galaxies have been a matter of debate. High resolution data from SPITZER, HERSCHEL 
and PLANCK allow us to probe the entire spectral energy distribution (SED) of morphologically separated components of the dust 
emission from nearby galaxies and allow a more detailed comparison between data and models. 

Aims. We wish to establish the physical origin of dust heating and emission based on radiation transfer models, which self-consistently 
connect the emission components from diffuse dust and the dust in massive star forming regions. 

Methods. NGC 4214 is a nearby dwarf galaxy with a large set of ancillary data, ranging from the ultraviolet (UV) to radio, including 
maps from SPITZER, HERSCHEL and detections from PLANCK. We mapped this galaxy with MAMBO at 1.2 mm at the IRAM 
30 m telescope. We extract separate dust emission components for the Hll regions (plus their associated PDRs on pc scales) and for 
the diffuse dust (on kpc scales). We analyse the full UV to FIR/submm SED of the galaxy using a radiation transfer model which 
self-consistently treats the dust emission from diffuse and SF complexes components, considering the illumination of diffuse dust both 
by the distributed stellar populations, and by escaping light from the Hll regions. While maintaining consistency with the framework 
of this model we additionally use a model that provides a detailed description of the dust emission from the Hll regions and their 
surrounding PDRs on pc scales. Due to the large amount of available data and previous studies for NGC 4214 very few free parameters 
remained in the model fitting process. 

Results. We achieve a satisfactory fit for the emission from Hll+PDR regions on pc scales, with the exception of the emission at 8 yum, 
which is undcrprcdictcd by the model. For the diffuse emission we achieve a good fit if we assume that about 30-70% of the emission 
escaping the Hll+PDR regions is able to leave the galaxy without passing through a diffuse ISM, which is not an unlikely scenario 
for a dwarf galaxy which has recently undergone a nuclear starburst. We determine a dust-to-gas mass ratio of 350-390 which is close 
to the expected value based on the metallicity. 

Key words, dust, extinction - Galaxies: individual: NGC 4214 - Galaxies: dwarf - Galaxies: ISM - Galaxies: star formation - 
submilUmeter: galaxies 



1. Introduction 

Interstellar dust is a component of the interstellar medium (ISM) 
that is present in all phases from dense molecular clouds to the 
warm ionised regions around massive stars. It plays an impor- 
tant role in star formation (SF) and in the overall energy budget 
of a galaxy. Measurements of the dust re-emission from galaxies 
is a powerful technique to quantify the SF. The details however, 
are not fully understood yet as the emission depends on sev- 
eral parameters as the geometry, the mixture between dust, stars 
and Hn regions and also on the dust properties in the interstellar 
medium which are expected to vary as a function of metallicity 
and are known to depend on the environment. 

The study of dust properties (e.g. mass, extinction coeffi- 
cient) has been an object of debate. Much progress has been 
achieved in recent years due to data from the satellites SPITZER, 
HERSCHEL and PLANCK, which allow for the first time to 
probe the entire wavelength range of the spectral energy distri- 
bution (SED) of the dust emission, from the mid-infrared (MIR) 
to the submillimeter (submm), for a large number of galaxies at 
the best resolution ever. 



The dust SED is a key observation for understanding the 
properties of dust. The dust SED of dwarf galaxies frequently 
shows differences to those of spiral galaxies. The main differ- 
ences are: i) The SED of dwarf galaxies shows a relatively low 
emission at 8 fim, most likely due to a relatively lower content 
of polycyclic aromatic hydrocarbons (PAHs) at low metallici- 
ties (e.g. Draine et al. 2007, Galliano et al. 2008, Engelbracht 
et al. 2008). ii) A submm "excess" has been found in the SED 
of many, mostly actively star-forming, low-metallicity galax- 
ies (Lisenfeld et al. 2002, Galliano et al. 2003, 2005, Bendo 
et al. 2006, Galametz et al. 2009, 2011, Israel et al. 2010, Bot 
et al. 2010, Dale et al. 2012, Planck Collaboration et al. 2011a). 
Several reasons have been suggested to explain this excess: the 
existence of a large amount of cold (< 10 K ) dust (Galliano et al. 
2003, 2005, Galametz et al. 2009, 2011), different dust grain 
properties ( Lisenfeld et al. 2002, Reach et al. 1995, Meny et al. 
2007, ), or magnetic nanograins, which can produce magnetic 
dipole radiation at microwave and submillimeter wavelengths 
(Draine & Hensley 2012). 

It is also noteworthy that one sees a variety of FIR colours 
for gas rich dwarf galaxies, with examples of warm FIR SEDs 
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(as found from IRAS colours - e.g. Melisse & Israel 1994), or 
of cold FIR SEDs (as revealed by ISO and Herschel measure- 
ments extending longwards of 100 fim - e.g. Popescu et al. 2002, 
Gross! et al. 2010). Such variety may not be too surprising given 
that the star formation in dwarf galaxies is likely fundamentally 
bursty in nature. If we interpret the FIR SED in terms of the com- 
bination of cooler distributed cirrus dust emission and warm dust 
emission from grains in SF regions, the latter opaque structures 
should be most prominent in the early stages of the evolution of 
a starburst leading to a warm FIR SED. over time the starburst 
fades, leaving the cirrus component and a cooler FIR SED. Other 
factors potentially contributing to the observed variety in FIR 
colours are differences in the contribution of very small grains 
(e.g. Galliano et al. 2003, 2005) and the possible presence of ex- 
tended cold dust outside the main star forming disk of the galaxy 
(Popescu et al. 2002). 

In order to interpret the dust SED of a galaxy and to under- 
stand the differences in the SEDs of dwarf galaxies, a physical 
model based on realistic dust properties and taking into account 
the heating and emission of dust immersed in a wide range of 
interstellar radiation fields (ISRFs) is needed. Ideally, radiation 
transport in a realistic geometry should be done, but this is of- 
ten difficult due to the complex geometry and large number of 
parameters. Models can generally be classified into three broad 
groups: (1) modified blackbody fits, which are too simple to de- 
scribe reality correctly but give a first idea of the dust temper- 
ature ranges, (2) semi-empirical models that try, in a simplified 
way, to describe dust immersed in a range of different ISRFs 
(e.g. Dale et al. 2001, Draine et al. 2007, Galametz et al. 2009, 
2011, da Cunha et al. 2008, Natale et al. 2010) and (3) mod- 
els that include full radiation transfer (e.g. Popescu et al. 2011 
for spiral galaxies, Siebenmorgen & Kriigel 2007 for star-burst 
galaxies, see also Silva et al. 1998, Popescu et al. 2000, Popescu 
et al. 2004, Misiriotis et al. 2001, Bianchi 2008, Baes et al. 2010, 
2011 and MacLachlan et al. 2011), which are the most precise 
description of a galaxy if all parameters, including the geometry, 
are known. 

The new IR facilities with their high angular resolution and 
sensitivity allow for the first time a detailed comparison between 
models and data. For nearby, and thus spatially resolved, galax- 
ies, the spatial variations of the dust SED can be probed and 
modelled. In particular, a good spatial resolution allows us to 
separate and treat the emission from dust heated by the UV ra- 
diation of massive stars in HII regions and their adjacent pho- 
todissociation regions (PDR), and the diffuse dust heated by the 
general ISRF. 

We choose NGC 4214 as a test case because of the 
large amount of ancillary data and its proximity (l"~14pc, 
D=2.9 Mpc; Maiz-Apellaniz et al. 2002), which allows us to ob- 
servationally separate the dust emission from the two prominent 
central SF complexes and from the diffuse dust component in 
the disk. This allows us to apply a radiation transfer treatment 
constrained by measurements of direct light in UV to NIR and 
dust/PAH re-radiated light in both the diffuse and SF complexes. 

We choose the radiation transfer model of Popescu et al. 
(2011) for spiral galaxies since this self-consistently treats the 
dust emission from the diffuse and SF complex components, 
considering the illumination of diffuse dust both by the smoothly 
distributed stellar populations, and by escaping fight from the 
HII regions in spiral disks. Although NGC 4214 is a dwarf 
galaxy, rather than a spiral galaxy, it nevertheless exhibits the 
principle geometrical features of the Popescu et al. (2011) 
model. In particular, NGC 4214 shows an exponential disk dis- 
tribution of the stellar light, the diffuse dust emission and the 



atomic gas emission, so is actually a good candidate to be fitted 
with a model for disk galaxies. While maintaining consistency 
with the framework of Popescu et al. (2011) model, we use the 
model of Groves et al. (2008) to provide a detailed description 
of the dust emission from the HII regions and their surrounding 
PDRs on scales of tens of pc. These PDRs mark the transition 
from the ionised medium to the dense molecular gas left over 
from the highly opaque cloud out of which the stars formed. The 
model of Groves et al. (2008) takes into account the dynamical 
evolution of the distance of the PDRs from the central ionis- 
ing star clusters due to the mechanical effect of stellar winds on 
the surrounding ISM. This model is therefore particularly suit- 
able to describe the two prominent SF complexes in the centre 
of NGC 4214, since these show shell structures imaged by HST 
surrounding the two main star clusters of the galaxy, which may 
be delineating wind blown bubbles around the HII regions. 

NGC 4214 is a Magellanic starbursting dwarf irregular 
galaxy (de Vaucouleurs et al. 1991a) which shows a large degree 
of structure, from HI holes and shells (Mclntyre 1998) typical of 
dwarf galaxies (Walter & Brinks 1999) to clear indications of a 
spiral pattern and a central bar. The molecular gas traced by the 
CO(l-O) line shows three well differentiated CO emitting com- 
plexes (Walter et al. 2001) related to the main SF complexes. 
NGC 4214 is a gas rich galaxy: the total mass of atomic gas is 
Mhi = 4.1 X 10^ Mq (Walter et al. 2008), and the molecular gas 
mass is Mm = 5.1 x 10* M© (Walter et al. 2001, obtained with 
a Galactic conversion factor). 

Karachentsev et al. (2004) reported for NGC 4214 a steUar 
mass of ~ 1.5 X 10^ Mq, similar to the mass found in the Large 
Magellanic Cloud (LMC). Ultraviolet (UV), optical and near- 
infrared (NIR) images of NGC 4214 show that the young stellar 
population is embedded in a smooth disk of old stars (see the 
large field of Fig. 1), which can account for a significant fraction 
(~ 75%) of the total steUar mass (Williams et al. 201 1). Despite 
this high fraction of old stars, NGC 4214 is a galaxy with an 
intense recent star formation activity, as shown by the two star 
forming complexes located in its centre and frequently referred 
in the literature as NGC 4214-NW and NGC 4214-SE. When re- 
solved, the two complexes show individual smaller knots of star 
formation (see the small field of Fig. 1). One of the most striking 
features in these complexes is the large shell structure in the NW 
region, where most of the gas in front of the central star cluster 
seems to have been removed by the action of stellar winds and 
supemovae (SNe) (Maiz-Apellaniz et al. 1998; MacKenty et al. 
2000). NGC 4214-SE is more compact and shows no clear ev- 
idence of a decoupling between the star clusters and the gas. 
The morphological differences for the two complexes were in- 
terpreted as an evolutionary trend by MacKenty et al. (2000). 
Using stellar synthesis models, Ubeda et al. (2007) determined 
the age, the mass, the radius and the extinction of the star clus- 
ters within the NW and SE complexes. They found an age of 
5 Myr for the star clusters in the NW region, whereas the age of 
the star clusters of the SE region ranges from 1.7 to 4.0 Myr. 

The intemal extinction distribution of the two complexes was 
studied by Maiz-Apellaniz et al. (1998). These authors found 
that the distribution of the dust clouds is correlated with the 
distribution of the ionised gas: the dust is located at the border 
of the star clusters in the NW complex, whereas for the SE re- 
gion the dust clouds seem to be co-spatial with the star clusters. 
The metallicities of the two complexes have been measured by 
Kobulnicky & Skillman (1996), who found values in the range 
of Z ~ 0.3 Zq with little dispersion, similar to the metalUcity of 
the Small Magellanic Cloud (SMC). 
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Thus, the large number and the broad range of previous stud- 
ies performed for this galaxy, as well as the significant amount 
of multiwavelength data, give us the opportunity of obtain an 
accurately picture of the stellar, dust and gas components in 
NGC 4214. 



2. The data 

A wide range of data exists in archives for NGC 4214, from 
the UV to the submm. The entire dust emission SED from 
the Mid-IR to the submm is covered by data from SPITZER, 
HERSCHEL and PLANCK. In addition, to extend the dust SED 
to the mm range we have mapped the galaxy at 1.2 mm with 
the IRAM 30 m telescope. In the following subsections we de- 
scribe the data and the reduction that we have performed prior 
to extracting the photometry. A selection of images at different 
wavelengths is shown in Fig. 2. 

2.1. GALEX 

NGC 4214 was observed with GALEX the 13th of January 
2010 as part of the program GI4-095 (Janice Lee, 2009). Far- 
ultraviolet (FUV) and near-ultraviolet (NUV) photometric maps 
were obtained. 

The GALEX FUV detector has an effective wavelength of 
1539 A, a bandwidth of 442 A and an image resolution of 4.2". 
For the NUV detector, the effective wavelength is 23 16 A, the 
bandwidth 1060 A and the image resolution 5.3". With both de- 
tectors the two main SF complexes NW and SE are well resolved 
and the field of view, with a radius of about 36', cover the entire 
disk of NGC 4214. We have retrieved the calibrated maps from 
the archive and subtracted the background using the correspond- 
ing maps from the archive. 

2.2. HST-WFC3 

Wide Field Camera 3 (WFC3) data were obtained for NGC 4214 
as part of the WFC Science Oversight Committee (SOC) Early 
Release Science (ERS) program (program ID 11360, P.I. Robert 
O'Connell). A total of 7 stellar and 7 nebular images were taken 
using UVIS and IR cameras. Both cameras onboard the HST 
have a field of view of ~2'and they cover the central star forming 
complexes NW and SE, allowing us to separate them into sev- 
eral smaller star forming regions. We obtained calibrated maps 
stellar filters F225W, F336W, F438W, F547M, F814W, Fl lOW 
and F160W from the HST Data Archive. 

2.3. SPITZER 

The SPITZER data for NGC 4214 used in this work are part of 
the Local Volume Legacy (LVL) survey (Dale et al. 2009). The 
LVL sample contains 258 galaxies within 1 1 Mpc, which have 
been mapped with both MIPS (3 bands) and IRAC (4 bands). 

For MIPS data the exposure times were 146.8, 83.8 and 
16.76 seconds at 24, 70, and 160 jum, respectively. The area 
mapped in the three bands, approximately 16'xl6', covers the 
whole disk of NGC 4214. The resolution of 6" for MIPS 24fj.m 
allows us to determine fluxes individually for the SF complexes 
NW and SE. This separation is not possible for MIPS 70 and 
160 jum due to a poorer (18" and 38", respectively) angular res- 
olution. We obtained the images from the archive and subtracted 
the background, which we determined as the mean value of an 
annulus located outside the disk of the galaxy. 



As for several other galaxies in the LVL sample, IRAC data 
of NGC 4214 was taken from previous SPITZER programs. In 
the case of NGC 4214, flie LVL makes use of the IRAC data 
from the "Mid-IR Hubble Atlas of Galaxies" (Fazio & Pahre 
2004). As part of this program, NGC 4214 was observed with 
IRAC at 3.6, 4.5, 5.8 and 8.0 /im in May 2004. In our study of 
the SF regions, we only use the IRAC 3.6/im band, which we 
assumed as pure stellar emission, and the IRAC S.Ofim band. 
In order to isolate the dust component from the IRAC 8.0 /im we 
subtracted the stellar emission using the IRAC 3.6 /im image and 
applying the formula provided in Helou et al. (2004). With an 
angular resolution of 1.7" and 2.0" for the IRAC 3.6 and 8.0 Aim 
bands, respectively, the SF complexes NW and SE are resolved 
into several smaller HII regions. 

The total mapped area by IRAC, approximately of 13' x6', 
covered the central part of the disk of NGC 4214, but not the en- 
tire extended disk. Therefore, the total emission at 3.6 and 8 jim 
that we obtain might be underestimated. However, judging from 
the IRAC 8.0 /im image in Fig. 2, we believe that the underesti- 
mate is not significant. 

The uncertainty in the flux calibration is better than 10% for 
aU IRAC (Reach et al. 2005, Fazio et al. 2004) and MIPS 24 j^m 
bands (Engelbracht et al. 2007), and is better flian 20% for MIPS 
70 and 160 //m bands (Rieke et al. 2004). We adopt these num- 
bers as the calibration uncertainties in our error estimates for the 
fluxes (see Sect. 3.5). 

2.4. HERSCHEL 

NGC 4214 was observed by HERSCHEL's PACS and SPIRE in- 
struments as part of the Dwarf Galaxy Survey programme (PI. S. 
Madden), a guaranteed time key program with the aim of map- 
ping the dust and gas in 51 nearby dwarf galaxies. 

SPIRE data were obtained the 26th of June 201 1 and cover a 
circular area of approximately 17'. We retrieved the HIPE Level 
2.0 images from the archive and performed further analysis for 
our purposes. Due to their good quality and absence of arte- 
facts, we only subtracted the background which we measured 
in the same annulus used for MIPS data. The resolution of 18" 
at SPIRE 250 fim allowed us to measure individual fluxes for the 
two star forming complexes. This distinction was not possible in 
the case of SPIRE 350 and 500 //m due to a poorer image resolu- 
tion (25" and 37", respectively). Swinyard et al. (2010) reported 
for SPIRE 250, 350 and 500 /xm an uncertainty in the flux cali- 
bration of 15%, which we adopted in our error calculation. 

PACS data were obtained the 27th and 28th of December 
2011 and cover an area of ~25'x25'. The data from the archive 
were processed to level 1 with HIPE 8.2.0. To complete the re- 
duction we used the Scanamorphos (Roussel 2012) map making 
technique. The background was subtracted from the final maps 
using the same annulus mentioned before. At the resolution of 
PACS at 70, 100 and 160 jum (5.2", 7.7" and 12", respectively), 
the two star forming complexes are resolved. The uncertainty in 
the flux calibration is of the order of 10%. 

2.5. PLANCK 

The "PLANCK Early Release Compact Source Catalogue" in- 
cludes fliree detections of NGC 4214 at 350, 550 and 850 /im, 
aU of them composed of two different observations. For the 
three detections the subtracted Cosmic Microwave Background 
(CMB) flux was less than 30% of the total flux measured orig- 
inally by PLANCK. The catalogue gives two estimates for the 
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Fig. 1. Combined SLOAN and HST image of NGC 4214. Large field: Image of NGC 4214 from a RGB combination of SLOAN 
r, g and u bands. Zoom-in: Image of the centre of NGC 4214 from a RGB combination of HST-WFC3 filters F657N (H,,, + cont), 
F502N ([OIII] + cont) and F336W (cont). The two circles are the apertures that we use for the aperture photometry of the regions 
NW and SE. 

total flux of a source: i) FLUX, obtained by aperture photom- 
etry within the nominal sky-averaged FWHM, which is 4.45', 
4.71' and 4.62' for 350, 550 and 850 /im, respectively, and ii) 
FLUXDET, obtained by their native detection algorithm (Planck 
Collaboration et al. 201 lb). We adopt here the first flux estimate, 
FLUX, which is recommended for sources that are point-like 
with respect to the PLANCK beam, as is the case for NGC 4214. 
We use the difference between the two estimates (which is larger 
than the nominal error of each flux) as an estimate of the errors. 

2.6. I RAM 30 m 

Several observations at 1200yL(m were made by our group 
between December 2009 and November 2010 at the IRAM 
30 m telescope on Pico Veleta (Spain), with the 37-channels 
bolometer array of the Max-Planck-Institut fiir Radioastronomie 
(MPIfR). The 37 pixels are located in a hexagonal structure with 
a beam-size of 10.8" and a pixel-to-pixel separation of about 
23". The observations were done on-the-fly and they were cal- 
ibrated by observations of the planet Mars and secondary cali- 
brators. Observations were reduced and combined into an image 
with Robert Zylka's MOPSIC' pipeline in a standard manner in- 
cluding baseline subtraction, spike and sky noise removal. We 
smoothed the combined image to a resolution of 18" in order 
to increase the signal-to-noise ratio. At this resolution, the two 
main star forming complexes are still resolved. Our observations 

' See http://www.iram.es/IRAMES/mainWik;i/CookbookMopsic 



were not sensitive enough to detect any extended emission out- 
side the two main SF complexes. We adopt 30% calibration un- 
certainty. 

2. 7. Thermal radio emission 

The thermal radio emission in an HII region is proportional to 
the production rate of Lyman continuum photons A^Lyc • Condon 
(1992) gives the following relation between the thermal ra- 
dio emission and the extinction-corrected Hfi line flux, F{lifi), 
which is proportional to A^Lyc^ 

1 / T f'^l V \-»' F(m 
mJy ~ 0.28ll04K/ \GHz/ lO-'^erg cm-^s-' ' 

where 5r is the thermal flux density at the frequency v for an 
HII region with a temperature T. 

We use the measurements of the Ha emission from 
MacKenty et al. (2000) to calculate the expected thermal radio 
emission from NW and SE. They presented an Ha narrowband 
map of the centre of NGC 4214 obtained with the HST Wide 
Field and Planetary Camera (WFPC2). They used the spectro- 
scopic values of the Balmer ratio Ha/HyS presented by Mafz- 
Apellaniz (2000) to correct for dust extinction and they reported 
intrinsic Ha fluxes of NW and SE. Assuming that the dust is 
well mixed with the gas, MacKenty et al. (2000) calculated 
extinction-corrected Ha fluxes of 12.88 x 10"'^ erg cm"^ s"' 
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for the NW complex and 5.75 x IQ-^^erg cm-^ for the SE 
complex. 

For the NW region, MacKenty et al. (2000) used a circu- 
lar aperture of radius 20.6", practically identical to the aperture 
used in this work (see Sec. 3.1). For the SE region they used a 
smaller aperture (9.96"). We therefore multiphed their flux by 
a factor 1.4, corresponding to the ratio of the uncorrected Ha 
fluxes measured from the same image in apertures of 9.96" and 
18", respectively. Using these extinction-corrected Ha fluxes, 
the intrinsic Ha I Hp ratio and the temperature of lO'^K reported 
by Kobulnicky & Skillman (1996), we obtain with Eq. 1 flux 
densities at 8.46 GHz (3.5 cm) of 13 and 8 mJy for NW and SE, 
respectively. 

The combined prediction for the thermal radio emission of 
NW and SE (21 mJy) can be roughly compared with the mea- 
surement from the VLA at 8.46 GHz of 24.2±4.8 mJy given 
by Kepley et al. (2011) for the total radio emission, which in- 
cludes synchrotron and thermal emission. From the uncorrected 
Ha emission, Kepley et al. (2011) estimated that the fraction of 
thermal emission is on average ~ 0.5 (see Fig. 9 in their paper). 
However, since these authors did not correct the Ha emission for 
dust attenuation, their estimation gives a lower limit of the ther- 
mal fraction. Our estimation of free-free emission (21 mJy) falls 
within the limits of thermal fractions of 0.5 (12.1 mJy) and 1.0 
(24.2 mJy) and is thus consistent with their data. We used their 
limits to roughly estimate an error of the thermal radio flux of 
~ 30%. 

3. Photometry 

We need to separate the dust emission from the two major SF 
complexes NW and SE and the diffuse disk, in order to estab- 
lish the individual SEDs of these morphological components. 
We therefore independently carried out aperture photometry for 
the regions SE and NW and for the total emission of the entire 
galaxy. We determined the diffuse dust emission as the differ- 
ence between the total emission and the sum of the emissions of 
SE and NW regions. 

Prior to extracting the photometry, we regridded all the im- 
ages to a common pixel size keeping the original resolution of 
each image. Then, corresponding aperture corrections were ap- 
pUed to the fluxes in each band. In this section we describe how 
we performed the aperture photometry, which corrections were 
appUed and how the errors were handled. 

3.1. Apertures 

We chose the apertures around the HII regions large enough to 
enclose the entire localized dust emission and simultaneously 
tried to include as little as possible emission from the diffuse 
component. Mainly based on the MIPS 24jum image, we chose 
apertures with radii of 21" and 18" for the NW and SE com- 
plexes (see Figures 1 and 2), respectively. The 24 jjm image was 
used as it is the filter which is expected to most faithfully trace 
the warm dust in thermal equilibrium with the strong radiation 
fields from the star clusters in the PDRs and HII regions associ- 
ated with the SF regions. We were able to measure the emission 
from each SF complex only for those wavelengths where both 
complexes were spatially resolved, which excluded the maps of 
MIPS 70 and 160 jum and SPIRE 350 and 500 ^um. The mea- 
sured, background- subtracted fluxes are listed in column 3 in 
Tables 1 and 2 for NW {F'^"^) and SE (F^'^''), respectively. 

One issue of particular relevance to the quantitative analy- 
sis presented in this paper is that of the subtraction of the local 



background underlying the NW and SE SF regions. Ideally, we 
need to cleanly separate the localised dust/PAH emission associ- 
ated with the HII regions and the local PDRs on scales of tens of 
parsecs (marking the interface between the tenuous gas ionised 
by the star cluster and the dense material of the parent molecu- 
lar clouds) from the underlying extended emission from diffuse 
dust distributed on the kpc scale of the disk. In practice it is very 
challenging to perform such a separation, not only because of 
blending effects due to the finite angular resolution (these are 
quantified through the aperture corrections described in Sect. 
3.2, below) but also because of the limited range of UV pho- 
tons which have escaped the SF regions and have passed into the 
diffuse dust layer, forming a halo around the SF regions. 

The extent and brightness of this diffuse halo emission is ex- 
pected to vary strongly with infrared wavelength in a way which 
is both difficult to measure directly and also uncertain to pre- 
dict theoretically on an a priori basis. Therefore, for simplicity, 
we have opted in this paper to make no attempt to subtract the 
underlying diffuse background when performing aperture pho- 
tometry on the SF regions. For future reference we simply note 
here that we may thereby be overestimating the flux densities of 
the NW and SE SF regions, most particularly in the PAH bands, 
but also to some extent around the peak of the SED in the FIR. 
The effect at 24 /im, however, is likely to be small, since, as al- 
ready noted, this band is primarily sensitive to warm dust heated 
by the intense radiation fields in the PDRs and HE regions close 
to the star clusters. 

For the measurement of the total emission of NGC 4214 we 
chose an aperture of 5' in radius (see Figure 2). We tested with 
growth curves at all wavelengths that this aperture enclosed the 
entire emission of the galaxy. For the PLANCK measurements, 
which include the total emission of the galaxy, the apertures 
could not be chosen. However, their values are very close (see 
Section 2.5) to ours, and our growth curves showed that at the 
radius of the PLANCK beam the emission was enclosed com- 
pletely. 

The total emission from the MAMBO 1.2 mm map is more 
uncertain due to the poor signal-to-noise ratio at the outer parts 
of the area covered by our observations, which is practically the 
same as the aperture of 5' used for the other images. In the mod- 
erately smoothed image (resolution 18", shown in Fig. 2) no dif- 
fuse emission can be seen. However, when we smooth the map to 
an angular resolution of 40" diffuse dust emission becomes vis- 
ible in the inner 160" whose structure corresponds very well to 
the SPIRE 500 //m map. Beyond this radius we are not sure that 
the structures are real. The growth curve of MAMBO keeps rais- 
ing until ~ 200" where it becomes approximately flat. We ob- 
tained integrated fluxes of 0.26±0.03 Jy and 0.35 ±0. 1 1 Jy for the 
apertures of 160" and 300", respectively. The lower value should 
represent a reliable lower limit for the total flux in NGC 4214. 
The total flux is more uncertain, however, given the flatness of 
the growth-curve beyond 200", we are therefore confident that 
the flux at 300" is a reasonable estimate. Due to these uncertain- 
ties, we did not include the MAMBO data point in our fitting 
procedure but we show it in the figures for comparison with the 
models. 

The measured, background- subtracted fluxes of the total 
emission are Usted in column 3 of Table 3. 

3.2. Aperture correction 

Because of the different resolutions of the images we applied 
aperture corrections to the NW and SE IR fluxes. We derived 
values of the aperture correction by adopting the IRAC 8 //m as 
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BAND 


Andinxi) 


^sub Hv) 


c 

^apert 


color 


^ NW v.*' J/ 


IRAC8 


7.872 


0.163 


1.000 


0.590 


0.276 ± 0.035 


M1PS24 


23.680 


0.733 


1.078 


0.986 


0.801 +0.103 


PACS70 


70.000 


7.612 


1.039 


0.982 


8.053 ±1.031 




1 nn nfvn 

iUU.UUU 




l.UJU 


yj,yoD 




PACS160 


160.000 


6.114 


1.075 


1.010 


6.507 ± 0.833 


SPIRE250 


250.000 


2.329 


1.110 


0.992 


2.605 ± 0.443 


MAMBO**' 


1200.000 


0.031 


1.067 


1.000 


0.033 ± 0.010 


St (0.5 cm) 


5000.000 


0.011 


1.000 


1.000 


0.011 ±0.003 



Table 1. Flux densities for region NW obtained for an aperture of 21" radius centred at RA=12hl5m39.6s Dec=+36°19'36."5 
(J2000) (see Figs 1 and 2). is the measured, background-subtracted flux and is the final flux, after applying aperture and 
colour corrections. The contribution from thermal radio emission was not subtracted. 



BAND 




F-" (Jy) 






FsE (Jy) 


1RAC8 


7.872 


0.105 


1.000 


0.590 


0.178 ±0.023 


MIPS24 


23.680 


0.602 


1.069 


0.986 


0.653 ± 0.084 


PACS70 


70.000 


4.659 


1.032 


0.982 


4.896 ± 0.627 


PACSlOO 


100.000 


5.416 


1.040 


0.985 


5.719 ±0.732 


PACS160 


160.000 


4.094 


1.057 


1.010 


4.285 ± 0.549 


SPIRE250 


250.000 


1.689 


1.078 


0.992 


1.835 ±0.312 


MAMBO<*) 


1200.000 


0.029 


1.021 


1.000 


0.029 ± 0.009 


St (0.5 cm) 


5000.000 


0.005 


1.400 


1.000 


0.007 ± 0.002 



Table 2. Flux densities for region SE obtained for an aperture of 18" radius centred at RA=12hl5m40.8s Dec=+36°19'05."5 
(J2000) (see Figs 1 and 2). is the measured, background- subtracted flux and Fse is the final flux, after applying aperture and 
colour corrections. The contribution from thermal radio emission was not subtracted. 



BAND 


/loOum) 


n?sK(Jy) 


C'color 


^DISK 


(Jy) 


Fdwf (Jy) 


IRAC8 


7.872 


0.678 


0.590 


1.149 ± 


0.115 


695-'° °" 


M1PS24 


23.680 


2.050 


0.986 


2.080 ± 


0.208 


625+" '°^ 


M1PS70 


71.440 


23.945 


0.893 


26.814 d 


: 5.363 


13 417+3 ''82 


PACS70 


70.000 


24.209 


0.982 


24.652 d 


: 2.466 


11.703!|:»^ 


PACSlOO 


100.000 


35.309 


0.985 


35.847 + 3.585 


21.128!!^^^ 


PACS160 


160.000 


34.097 


1.010 


33.759 ± 3.376 


22.967!J^«« 


MIPS 160 


155.900 


39.016 


0.971 


38.574 ± 7.715 




SP1RE250 


250.000 


18.888 


0.992 


19.033 d 


: 2.855 


14 594+2-20* 

2g4j 


PLANCK350 


350.000 


9.546 


0.986 


9.682 ± 


0.988 


7 Q^^+a^97 
' -0.870 


SPIRE350 


363.000 


10.059 


0.999 


10.068 ± 1.510 


8.453+|J8| 


SPIRE500 


517.000 


4.497 


1.025 


4.387 ± 


0.659 


3.848!«J^» 


PLANCK550 


550.000 


3.635 


0.921 


3.947 ± 


1.318 


3.510!!;™ 


PLANCK850 


850.000 


0.946 


0.887 


1.066 ± 


0.264 




MAMBO<*' 


1200.000 


0.350 


1.000 


0.350 ± 


0.106 


0.288!o;i 



Table 3. Flux densities of NGC 4214 for an aperture of 5' radius centred at RA=12hl5m39.1s Dec=-H36°19'34."5 (J2000) (see 
Fig. 2). -pQisf^ is the measured, background- subtracted flux, Fdisk is the final flux of the entire galaxy, after applying colour 
correction, and Fdiff is the flux of the difliise emission obtained as Fuisk - (^nw + ^se)- The lower error of Fdiff includes 
the contribution of smaller and less intense SF regions in the disk (see Sect. 6). The contribution from thermal radio emission 
was not subtracted. 
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the high-resolution model of the flux distribution. We first mea- 
sured the flux within the NW and SE apertures in the IRAC 8 /xm 
image. Then we convolved the IRAC 8jum image with the 
point-spread-function (PSF) of MIPS 24, PACS 70, PACS 100, 
PACS 160, SPIRE 250 and MAMBO, and measured again the 
flux within the same apertures. The ratio between the flux mea- 
sured in the original and convolved images gives the aperture 
correction for each band. The values are listed in column 4 of 
Tables 1 and 2. 

3.3. Colour corrections 

In order to directly compare the monochromatic fluxes with the 
models of Popescu et al. (2011) and Groves et al. (2008), we 
appUed colour corrections to flie IRAC, MIPS, PACS, SPIRE 
and PLANCK fluxes. We followed the procedures described in 
the corresponding manuals^, where colour correction factors for 
diff"erent source distributions are listed. 

We chose the colour correction factors for the spectrum that 
most closely resembled the model spectrum at a given wave- 
length. For the MIPS 24 jim band we adopted a blackbody with 
a temperature of 70 K, for the MIPS 70 //m, MIPS 160 fim and 
PACS bands a blackbody with a temperature of 50 K, for SPIRE 
bands and PLANCK 350 //m we chose a power-law spectrum 
with a = -1-2.0, for PLANCK 550 fim a power-law spectrum 
with a = +2.5, and for PLANCK 850 jum a power-law spectrum 
with a = -1-3.0 (where the exponent a is defined as Fy oc y"). 
For the particular case of IRAC 8 jum, we adopted the colour 
correction for a PAH-dominated spectrum. At 8 /xm we addi- 
tionally compared model and data by generating synthetic inte- 
grated fluxes. For this we (i) multiplied the model spectrum with 
the 8 fim filter profile and integrated it over the 8 fim filter band 
and (ii) multiphed the observed 8 jum data, adopting the black 
body of T = 10'' K profile on which the flux definition of IRAC 
is based, with the 8 fim profile and integrated it over the 8 /im 
filter band. The comparison of both values allowed an alternative 
estimate of the goodness of the fit, based entirely on the data and 
the model spectrum. We found that this comparison produced 
within 10% the same ratios as the comparison between colour 
corrected data point and model flux. 

The colour corrections are listed in Tables 1, 2 and 3. 

3.4. Line contamination 

Most of the FTR/submm broad-band filters used for this work 
are contaminated by molecular or atomic emission lines. In the 
cases where data were available, we estimated the line contam- 
ination by comparing the luminosity of the emission line with 
the total luminosity measured in the filter, obtained by integrat- 
ing the product of the source spectrum and the spectral response 
of the filter over the frequency range of the filter band. 

Emission lines from different CO rotational transitions fall 
into SPIRE, PLANCK and MAMBO filters. We used the CO(l- 
0) map presented by Walter et al. (2001) as a reference. We de- 
termined a velocity integrated flux of 9.3 Jy km s"^ in the NW 
region and 12.8 Jy km s ' in the SE region. The MAMBO fil- 
ter band, with a bandwidth of about 80 GHz, is affected by the 
CO(2-l) fine. We estimate the flux of flie CO(2-l) fine from 
the measured CO(l-O) value assuming that the hne intensity ra- 



^ see the IRAC Instrument Handbook, MIPS Instrument Handbook, 
PACS Observer's Manual, SPIRE Observers' Manual and Planck 
Collaboration (2011) for IRAC, MIPS, SPIRE and PLANCK, respec- 
tively. 



tio it is given by the value for optically thick, thermalised CO, 
^co(i-0)/^co(2-i) = 1- The contribution of the CO(2-l) line to 
the total measured flux at 1.2 mm of the NW and SE complex is 
then 0.87% and 1.27%, respectively. For other bands, the con- 
tamination by higher CO transitions is more than an one order 
of magnitude smaller. 

PACS spectroscopic maps of atomic FIR lines ([OI] 63, 
[OIII] 88, [Nil] 122, [OI] 146, [CII] 158 and [Nil] 205 yum) 
of NGC 4214 were presented by Cormier et al. (2010). With a 
field of view of 1.6' x 1.6', only the central part of the galaxy is 
observed. However, this region completely covers the main SF 
complexes NW and SE where the major part of the line emission 
originates. Therefore, we do not expect to severely underesti- 
mate the total emission of these lines. [OI] 63 and [OIII] 88jum 
fall into the bandpass of the MIPS 70 and PACS 70 filter. The lu- 
minosities reported by Cormier et al. (2010) are 0.89 x 10*" Lq 
for [OI] 63 jj.m and 1.99 x 10*^ Lq for [OIII] 88 //m. The sum 
of these values corresponds to 1.11% (1.33%) of the total lumi- 
nosity measured in the PACS 70 (MIPS 70) band. [OIII] 88jum 
fall into the bandpass of the PACS 100 filter. Its luminosity 
corresponds to 1.73% of the total luminosity measured in the 
PACS 100 band. [OI] 146 and [CII] 158 //m fafl into the bandpass 
of the MIPS 160 and PACS 160 filter. The reported luminosities 
for[OI] 146and[CII]//mUnesare0.05xl06Lo and2xl0^Lo , 
respectively, which together represents 2.88% (4.26%) of the to- 
tal luminosity measured in the PACS 160 (MIPS 160) band. 

Due to the small contribution from the atomic and molecu- 
lar lines to the FIR/submm bands, we did not apply any decon- 
tamination to our measurements except for the PACS 160 and 
MIPS 160 bands. These have been included in Tables 1, 2 and 3. 

3.5. Error handling 

In our error analysis we took into account three types of error: (i) 
calibration, Acai, for which we adopted the values mentioned in 
Sect. 2, (ii) measurement error due to background fluctuations, 
Aback, and (iii) an estimate of the error due to the uncertainty in 
the aperture size of regions SE and NW, Aapert (this error was not 
relevant for the total emission). We neglected the uncertainties 
introduced by the colour and aperture corrections. 

The error due to the background fluctuations was calculated 
assuming that each pixel within the aperture has an error given 
by the standard deviation of the background noise, o"back- In ad- 
dition, we had to take into account the error of the background 
which was subtracted within an aperture of A^apen pixels. This 
error is o'backA'apert/ VMjack, where A'back is the number of pixels 
used to compute the level of background. This gives (see also 
Dale et al. 2012) a total error for the background subtracted flux 
of: 



Aback - cr back A/A^apert + T^^ (2) 



We estimated the error Aapert by changing the aperture size 
for the regions NW and SE by ±1". We obtained changes in the 
integrated fluxes between 4 and 8% for the difl'erent bands. We 
adopted conservatively an error of Aapen - 8% for all bands. 

The final error for the flux is the quadratic sum of Acai, Aback, 
and Aaper- This is listed in Tables 1, 2 and 3 for the fluxes in each 
band. The dominant error sources were Acai and Aaper, whereas 
Aback was found to be neghgible for most wavelengths. 
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4. Models for the dust emission 

We analyse the full UV to FIR/submm SED of the different emis- 
sion components of the galaxy using the radiation transfer model 
of Popescu et al. (2011), which self-consistently treats the dust 
emission from diffuse and SF complexes components, consider- 
ing the illumination of diffuse dust both by the distributed stellar 
populations, and by escaping light from the HII regions. While 
maintaining consistency with the framework of Popescu et al. 
(201 1) model, we use the model of Groves et al. (2008) to pro- 
vide a detailed description of the dust emission from the two 
central SF complexes NW and SE. In the following sections we 
present a brief description of the physics and parameters of these 
models. 



4.1. The model of Popescu et al. (2011) 

Popescu et al. (201 1) presented a self-consistent model based on 
full radiative transfer calculations of the propagation of starlight 
in disk galaxies. To approximate the large scale geometry of the 
galaxy (see Fig. 1 in Popescu et al. 2011), they used two sepa- 
rate components: (i) an old component consisting of an old stel- 
lar disk, an old stellar bulge and a dust disk, and (ii) a young 
component consisting of a young stellar disk and a dust disk. 
The young component is introduced to mimic the more complex 
distribution of young stars and diffuse dust associated with the 
spiral arms. Apart from the diffuse component, the model in- 
cludes a clumpy component, consisting of the parent molecular 
clouds of massive stars. The input parameters of this model are: 

- The total central face-on B-band opacity, r^, which is the 
sum of the central face-on B-band opacities of the yoimg and 
old dust disks. 

- The star formation rate, SFR. 

- The dumpiness factor F, which can be physically identi- 
fied with the luminosity-weighted mean fraction of direc- 
tions from the massive stars, averaged over the lifetime of the 
stars, which intersects the birth-cloud. The dumpiness factor 
F is linked to the fraction of photons that escapes (/esc) from 
the SF regions into the diffuse medium as F = l-/esc- 

- The normalized luminosity of the old stellar disk, old. 

- The bulge-to-disk ratio, B/D, which determines the bulge 
contribution to the old stellar radiation field. 

- The radial scale-length of the old stellar disk, h,, which de- 
fines the size of the galaxy. All other spatial scales in the 
galaxy (scale-length of the young stellar population, scale- 
length of the dust and vertical scale-heights for the different 
components) have a constant ratio with hg (see Table E. 1 in 
Popescu et al. 2011). 

- The inclination angle, i. 

From the primary parameters SFR and F, Popescu et al. 
(201 1) define the star formation rate powering the diffuse emis- 
sion, SFR', as it follows (Eq. 45 in Popescu et al. 201 1): 



SFR' =SFRxil- F) 



(3) 



The Ubrary of diffuse SEDs of Popescu et al. (201 1) contains 
results for a four-dimensional parameter space spanned by t^, 
SFR', old and B/D. The diffuse component is calculated as an 
extrinsic quantity corresponding to a reference size (correspond- 
ing to a reference scale-length). In order to scale the intensity 
of the ISRF heating the diffuse dust in a galaxy, the parameters 



S FR' and old must be scaled to the reference size by comparing 
the scale-lengths (Eq. D.3 in Popescu et al. 201 1): 



SFR 



model 



SFR' X 



jhfiB) 
hs{B) 



old 



model 



= old X 



hs{B) 



(4) 



(5) 



where li'^iB) = 5670 pc and h,{B) are the reference B-band 
scale-length and the B-band scale-length of the galaxy under 
study, respectively. We note here that S FR'""'^"' and old'""'^"' are 
only internal parameters that allow us to interface with the li- 
brary of models. Because of this all the results in this paper are 
presented only in terms of the S FR and old, the real parameters 
of the galaxy under study. 

An additional scaling is required to set the flux levels of the 
SEDs from the library to the observed SED of our galaxy. Thus, 
the SED that represents our galaxy, Ff, is detennined as (Eq. D.2 
in Popescu et al. 201 1): 



hAB) 
yhf{B) 



. pd,model^^ijj^ f ^ pj^model^ old'""''''). (6) 



We would like to stress that in these models the absolute flux 
level of the predicted dust SED is fixed by the input parameters. 

The models of Popescu et al. (2011) have been designed 
to provide the formalism for fitting the integrated emission of 
galaxies, since for most galaxies resolved information on the 
dust emission is not available. In particular dust emission SEDs 
of individual star-forming complexes cannot be derived in most 
cases. Consequently, the model of Popescu et al. (201 1) provides 
an average template SED for the dust emission of the ensemble 
of the HII regions. This template has been empirically calibrated 
on data of a representative sample of prominent star forming 
complexes in our Galaxy, as fitted using the model of Groves 
et al. (2008). By contrast, in tiiis study of NGC 4214, we have 
resolved information on the central star forming complexes, and 
can therefore replace the average template with a detailed mod- 
elling of these complexes using directly the model of Groves 
et al. (2008), while still retaining the general framework of the 
model of Popescu et al. (2011), when calculating the attenua- 
tion of stellar light in the clumpy component. This allows us to 
reach a self-consistent treatment of the diffuse and clumpy com- 
ponent, mediated by the fact that the two models used here have 
a common parameter, the F factor. 

4.2. The model of Groves et al. (2008) 

The model of Groves et al. (2008) describes the luminosity evo- 
lution of a star cluster of mass Ma, and incorporates the expan- 
sion of the HII region and PDR due to the mechanical energy 
input of stars and SNe. The dust emission from the HII region 
and the surrounding PDR is calculated from radiation transfer. 
The main input parameters of the model are: 

- The metallicity of the star cluster, Z, in units of the solar 
metallicity Zq (Asplund et al. 2005). The metallicities used 
in the model are restricted to the SMC and LMC metallic- 
ities. This limitation is introduced by the stellar population 
model (Smith et al. 2002). 

- The age of the star cluster, t. 
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- The ambient pressure, log{po/k), which controls the rate of 
the expansion of the HII region and PDR due to the mechan- 
ical energy input of stars and SNe, and the size. Therefore 
\og(po/k) determines the dust temperature for a given star 
cluster mass and age. 

- The compactness parameter, log(C), which parametrizes the 
heating capacity of the star cluster and depends on Ma and 
log(po/^)- 

- The covering factor, /cov, which represents the fraction of 
the surface of the Hll region covered by the PDR. This is the 
same parameter as the F factor in the model of Popescu et al. 
(2011). 

- The hydrogen column density of the PDR, Nl 



PDR 



At a given age, the model self-consistently calculates the lu- 
minosity distribution of the star cluster and the radius of the inner 
edge of the PDR and computes the emergent SED for the dust 
in the HII region without considering the PDR, F^", and for the 
dust from the Hll region which is completely covered by a PDR, 
pHii+PDR rpjjg iQ^\ emergent SED from the star forming region 

which is partially covered by the PDR, F^^, is then given by: 



f'f = (l-/-v)xFf^ + /c 



V F' 

cov X 



HII+PDR 



(7) 



5. Constraints on thie input parameters 



The large amount of ancillary data covering a wide wavelength 
range, as well as the result of previous studies from the litera- 
ture allow us to determine or at least constrain most of the input 
parameters for both models. 

5.1. Input parameters for the model of Groves et al. 2008 

The two main SF complexes NW and SE are formed by an en- 
semble of smaller HII regions. Ideally, aU these smaller knots 
of SF should be modelled separately. However, IR/submm ob- 
servations lack the angular resolution to perform such a detailed 
analysis. Due to this hmitation, for both complexes we adopted 
parameter ranges which are wide enough to enclose the values 
of their individual knots. With the only exception of N™^, for 
which we adopted the value 10^^ cm"^, typical of molecular 
clouds in our own galaxy, the input parameters of the model of 
Groves et al. (2008) were observationally constrained: 

(i) Metallicity: The metaUicities of both SF regions have 
been measured by Kobulnicky & Skillman (1996), who found 
values of 12+log(0/H) = 8.17± 0.02 for the NW and 8.27+ 0.02 
for the SE region. With the Solar abundance 12+log(0/H) - 8.66 
(Asplund et al. 2005), as used in Groves et al. (2008), this gives 
Z = 0.32 Zo for NW and Z = 0.41 Zg for SE. We used a metal- 
licity of Z=0.4 Zg and also tested Z=0.2 Zq, which yielded a 
similar result (no template for Z=0.3 Zq is available in Groves 
et al. (2008)). 

(ii) Age: The age of the star clusters of NGC 4214 have 
been extensively studied by different authors. Leitherer et al. 
(1996) presented HST Faint Object Spectrograph (EOS) ultra- 
violet spectra covering the main star cluster of the NW complex. 
Using spectral synthesis modelling, Leitherer et al. (1996) re- 
ported for this stellar cluster an age of 4-5 Myr. Maiz-Apellaniz 
et al. (1998) obtained optical long-sht spectra of NGC 4214 
with the ISIS spectrograph of the William Herschel Telescope. 
Their bidimensional spectra covered both SF complexes. From 
the comparison of different observationally determined param- 
eters (equivalent width of Hfi, Wolf-Rayet population, effective 



temperature and UV absorption lines) to synthesis models, they 
determined ages of 3±1 Myr for both regions. MacKenty et al. 
(2000) mapped the central part of NGC 4214 with the HST 
Wide Field Planetary Camera 2 (WFPC2). From the Ha equiv- 
alent width, these authors determined an average age of 3.0- 
4.0 Myr and 2.5-3.0 Myr for the NW and SE regions, respec- 
tively. More recently, Ubeda et al. (2007) also used the WFPC2 
data set to determine the ages of the star clusters by a likelihood- 
maximization technique from the photometric colours. They de- 
termined an age of 5 Myr and of 2-4 Myr for the star clusters 
within the NW and the SE complex, respectively. Based on the 
dispersion found by these studies, we adopted ages of the NW 
and SE complexes of 3-5 Myr and 2-4 Myr, respectively. 

(Hi) Ambient pressure and Compactness: log(po/k) and 
log(C) were determined by comparing the expected and the ob- 
served radii of the individual HII regions as a function of the 
age. For a cluster of a given age it is possible to find different 
combinations of log(pf)/k) and log(C) that provide the observed 
radius of the expansion bubble. However, this degeneracy can be 
avoided when the mass of the star cluster is known, since both 
parameters log(po/^) and log(C) are related by the equation: 



log 



Me, 

Mo 



xlog(C) 



2^, / Po/k 
3 \cm 



(8) 



From the values of the masses reported by Ubeda et al. 
(2007) we obtained that for the NW complex the value of 
\ogipo/k) ranges from 7.0 to 8.0 and log(C) ranges from 5.5 to 
6.5 depending on the HII region considered. For the case of the 
SE complex, we found that log(po/^) ranges from 6.0 to 7.0 and 
that log(C) ranges from 4.5 to 5.5. We used these intervals to 
constrain both parameters. 

(iv) Covering factor: We assumed that the PDR consists of 
optically thick, homogeneously distributed clouds that surround 

the star cluster, leaving a fraction uncovered so that the light 
can escape unattenuated from these "holes". To a good approx- 
imation, the intrinsic luminosity of the central star cluster is the 
sum of the observed luminosities of the stars, Lstar. and the lumi- 
nosity re-emitted by the dust, Z-dust- The covering factor is then 
/cov = idust/(idust + ismr). We obtained Ldust by integrating the 
best-fit model template (see Sect. 6) from 3 fim to 1 .5 mm. To 
calculate Ls,ar we integrated the de-attenuated fluxes measured 
in our apertures for the SF regions from GALEX FUV to IRAC 
3 fim. The exact value of the opacity in front of the SF regions 
due to the diffuse dust layer is not known. For this reason, we 
carried out two different estimates, the first based on the lowest 
and the second on the highest realistic opacity, i) We adopted the 
foreground opacity of Ty = 0.35 (Ubeda et al. 2007) measured 
locally in front of the main stellar cluster of the NW region and 
assumed that this value is representative for both SF regions. 
The stars in the NW region seem to have evacuated most of the 
surrounding material associated with the HII region so that this 
value is expected to be lower than the average disk opacity in 
this area. We obtained /^w = 0.45 and /fj^ = 0.65 for the NW 
and SE regions, respectively, ii) We used the value of = 2.0 
which is the best fit value derived from our modelling of the 
diffuse emission (see Sect. 6). Assuming that the regions are in 
mid-plane and using the inclination angle i = 44° (Walter et al. 
2008), the opacity in front of the SF complexes is calculated as 
0.5 X T^/ cos(i) = 1.4. In this case, we obtained = 0.20 and 
fcm - 0.30. Thus, we obtain estimates for the covering factor of 
/^* = 0.20 - 0.45 and f^^^ = 0.30 - 0.65 for both regions. 
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Fig. 3. B-band scale-length of the stellar disk of NGC 4214. 
The fit (solid line) was achieved using the isophotes in the ra- 
dial range marked by the grey-shaded area. The uncertainty in 
the slope is represented by the dark grey area and the dark grey 
dashed lines. 

5.2. Input parameters for the model of Popescu et al. (2011) 

In the case of the model of Popescu et al. (201 1) the parameters 
that we could determine observationally are: 

(i) Bulge-to-disk ratio: Since no bulge is visible in the 
late-type galaxy NGC 4214 we set B/D = 0. 

(ii) Inclination angle: We fix i = 44° following Walter et al. 
(2008). 

(Hi) Scale-length: We determined the radial stellar scale- 
length in the B-band from a SLOAN g-band image. First, we 
removed the background of the image by calculating the average 
value of a wide set of small circular apertures placed strategically 
outside the disk of NGC 42 14. Then, we removed the contamina- 
tion from foreground stars and Hll regions, including the central 
regions SE and NW. We used the IRAF task ellipse (STSDAS 
package) with steps of 20" to determine the isophotes. For all 
the isophotes beyond 100" from the centre, the PA and the el- 
lipticity were fixed. Finally, we fitted an exponential function to 
the surface brightness profile determined with the mean values 
in each isophotal ellipse. We considered the ellipses whose mean 
values and the corresponding uncertainties were above the back- 
ground level. 

We found that it is not possible to fit the profile for the whole 
disk with a single exponential function. For this reason we 
decided to calculate separately the scale-length in the inner and 
the outer part of the disk (see Fig. 3). We visually determined 
the change of slope to take place at a radius ~ 90". The 
profile of the inner part of the disk (excluding the HII regions 
NW and SE) can be well fitted with a scale-length of ~450 
pc. For the outer part we determined a scale-length of SVS^jjj pc. 

(iv) dumpiness: As a first estimate of the dumpiness factor 
F we used the mean value from the upper and the lower limits 
obtained for the two SF regions, i. e., F - 0.40 ± 0.20. This 
parameter is only used for the modelling of the total integrated 



UV/optical/FIR emission described in Sect. 6.3. In Sect. 6.2 we 
will present an additional global estimate of F. 

(v) old: We derived this parameter as old = 0.009 by assum- 
ing that all the luminosity in the J, H and K bands (see Table 4) 
corresponds to the old stellar population. We integrated the lu- 
minosity of NGC 4214 in the J, H and K bands and derived old 
as the ratio of this luminosity and the integrated luminosity of 
the model galaxy (Table E2 in Popescu et al. 201 1). 



BAND 




^DISK 


(Jy) 


Reference 


FUV 


0.154 


0.072 ± 


0.007 


1 


NUV 


0.232 


0.091 ± 


0.009 


1 


B 


0.445 


0.354 ± 


0.056 


2 


V 


0.551 


0.446 ± 


0.064 


2 


J 


1.220 


0.520 + 


0.014 


3 


H 


1.630 


0.614 ± 


0.023 


3 


K 


2.200 


0.458 ± 


0.022 


3 



Table 4. Total stellar flux densities of NGC 4214. (1) This work, 
(2) de Vaucouleurs et al. (1991b), (3) Jarrett et al. (2003). 



6. Results 

In a first step we separately fitted the emission from the SF re- 
gions and from diffuse dust. In a second step we self-consistently 
combined the emission from the SF regions with that from the 
diffuse dust to fit the total emission of the galaxy. 

The emission from the NW and SE regions have been di- 
rectly measured from the maps as described in Section 3. The 
diffuse emission has been calculated as the difference between 
the total emission and the sum of the emission from the two SF 
complexes. For the data points where no direct measurement was 
available, the best-fit model value (see below) was taken. The 
values for the diffuse dust emission are listed in Table 3. 

In order to take into account the emission from smaller and 
less intense SF regions in the disk of NGC 4214, we estimated 
their contribution from the Ha emission of the ten brightest sec- 
ondary SF regions which corresponds to 18% of the Ha emis- 
sion of SE+NW. We assume that the shape of the dust SED from 
these smaller SF regions is the same as the sum of the SEDs of 
SE+NW, and we include their contribution in the calculation of 
the lower error range of the fluxes of the diffuse dust component 
in each band (see Table 3). The smaller HII regions only have a 
noticeable effect on the emission at 24 fim and 70 fim. 

6.1. Best fits for the SF regions 

The best-fit models for the NW and SE complexes are shown in 
Figures 4 and 5, respectively. 

We obtained a good fit for all our observed photometric 
points longwards of lOjum. For a metallicity of Z=0.4 Zq the 
best-fit models of the two complexes correspond to the param- 
eters f= 4.0Myr, log(C)= 5.0, log(po/k)^ 8 and /cov= 0.30 
for both regions. With the exception of the log(C) of the NW 
region and log(p()/k) of the SE region, all the parameters fall 
within our parameter ranges constrained from the observations 
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Fig. 4. Best-fit models to the NW complex for Z=0.4 Zq (solid line) and Z=0.2 Zq (dashed line). The values of the reduced 
(neglecting the IRAC 8 fim point) are 0.25 and 0.75 for Z-0.4 and 0.2 Zq, respectively. 



10,000 - 



1.000 r 



<^ 0,100 



0,010 - 



0,001 r 




X (ttm) 



Fig. 5. Best-fit model to the SE complex for Z=0.4 Zq (solid line) and Z-Q.2 Zq (dashed line). The values of the reduced 
(neglecting the IRAC 8jum point) are 1.17 and 0.89 for Z=0.4 and 0.2 Zq, respectively. 
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(see Sect. 5.1). For a metallicity of Z-0.2 Zq the best fit param- 
eters are in the case of the NW region t= 5.0 Myr, log(C)= 5.0, 
log(po/k)= 1 and /cov= 0.30 and in the case of the SE region 
f= 3.5 Myr, log(C)= 4.5, log(/7()//t)= 8 and 0.60. The value 
for log(C) in NW is slightly lower, and logipo/k) in SE slightly 
higher, but the rest of the parameters are within the observed 
ranges. 

We note that the IRAC 8 //m data point was excluded from 
our fitting procedure. The reason is that both models consid- 
erably underestimate the emission at 8 ^m. Specifically, from 
our best-fit models for Z=0.2 Zq we obtained that the observed 
fluxes are larger than the model values by a factor of 3.8 and 4.0 
for NW and SE, respectively. In the case of the best-fit models 
of Z=0.4 Zq the discrepancy decreases to a factor of 3.1 and 
2.9 for NW and SE, respectively. A complete discussion of the 
IRAC 8 jum discrepancy will be given in Sect. 7. 

6.2. Best fit for the diffuse emission 

We searched the library of the diffuse dust SEDs of Popescu et al. 
(201 1) for the best fit to the data in the MIR/submm range, leav- 
ing and S FR' as free parameters and keeping old fixed to the 
value 0.009. Fig. 6 shows the best fit obtained for our measured 
value hs = 873 pc (solid red fine), where we determined as best- 
fit values = 2 and SFR' = 0.059 Moyr ^. Neither flie data 
point at 8 jim micron nor at 1.2 mm were taken into account in 
the fitting procedure, the latter due to the observational hmi- 
tations discussed in Sect. 3.1. 

The model fits the data points in general very well. The 
largest discrepancy occurs at 8 /urn where the model overpredicts 
the observations by a factor of 1 .9. We will discuss this further in 
Sect. 7. At 160 fim there is a discrepancy with the PACS value 
which is overpredicted by 34% by the model, but the model pre- 
diction agrees within the errors with the MIPS data point at the 
same wavelength. The model overpredicts the SPIRE 250 fim 
point by 24% and underpredicts the 850 fim point by 34%. For 
the other data points the model predictions agree within the error 
bars. 

It is instructive to compare this value of SFR' derived from 
the fit to the diffuse component of the dust emission, with the 
value for S FR' derived from our measurement of the spatially 
integrated flux density in the UV-to-blue band and which we will 
call, for clarity, S FR'uy in the following. Applying the attenua- 
tion corrections due to diffuse dust, as tabulated in Popescu et al. 
(2011) for Tg = 2 and an inclination of 44° we derive a value 
of Luv = 3.70 X 10^^ W. We then determine SFR'^y/ follow- 
ing Eq. 17 in Popescu et al. (2011) as SFR'uy = Luv/L^"""^^, 
where L,-'^^^""fjy is the normalization factor We obtain S FR'^y - 
0.165 M0yr"'. This is a factor of 2.8 times higher than the 
value of S FR' derived from the fitting of the diff'use dust emis- 
sion. The most straightforward explanation for this discrepancy 
would be that the true value for the disk scale length hs is higher 
than the adopted of 873 pc. This would favour a disk with lower 
dust surface densities, leading to a higher value of SFR' needed 
to account for the observed amplitude and colour of the dif- 
fuse dust emission. If we adopt for the scale-length the maxi- 
mum allowed by the B-band images, hs = 1045 pc, we obtain 
a best fit to the difi\ise dust emission SED of = 1.2 and 
SFR' = 0.088 Moyr"^ At the same time, this reduction in Tg 
means that the deattenuated, integrated UV luminosity, Luv, is 
now lower, yielding S FR'^y - 0.147 M0yr"' which is only a 
factor of 1.67 greater than the value of SFR' derived from the 



diffuse dust emission. A complete discussion of this discrepancy 
will be given in Sect. 7. 

As a further consistency check we can use the value of S FR' 
derived from the diff'use dust emission to make a different, global 
estimate of the fraction of the UV radiation escaping from the 
SF regions, /esc, and the corresponding dumpiness factor, F, 
and compare it with the corresponding value already derived in 
Sect. 5 from analysis of the direct and dust-re-radiated UV light 
seen from the spatially resolved SF regions. For this, we have 
to assume that the effective star formation rate powering the dif- 
fuse dust, SFR', is due to radiation escaping from the regions 
NW and SE. The total effective SFR produced by NW+SE is 
then S FR'+S FRioc where SFRioc is the effective star-formation 
rate needed to power the dust emission from the SF regions. We 
calculate 5F/?ioc by integrating the dust SED from the SF re- 
gions over the entire wavelength range, obtaining Ldust, and as- 
suming that this luminosity is equivalent to the same amount of 
UV radiation absorbed by the dust locally. We then determine 
S FR\oc in an equivalent way as described for SFR'^j^ above, as 
SFRi,, = Li^jUZlv We obtain Lj^^t = 9.06 x lO^^ w, which 
gives SFR^^^ = 0.040 Moyr"'. With SFR' = 0.059 Moyr"', 
derived from our best fit for the diffuse emission, we then obtain 
/esc= SFR' I {SFR' +SFRioc) = 60% of the UV radiation of the 
SF region is required to escape from the SF complexes in order 
to heat the diffuse dust. This escape fraction corresponds to a 
/cov = F = I - 0.60 = 0.40, which is the same value derived in 
Sect. 5. We thus adopt in the following F = 0.40 ± 0.20 (where 
the error comes from the range of /cov derived in Sect. 5.1). 
From the best-fit value of SFR' and F = 0.40 we can esti- 
mate the value of the total star formation rate of NGC 4214 as 
SFR = SFR'/il -F)^ 0.059/(1 - 0.40) = 0.098 Mqyt-K 

6.3. Best fit for the total emission 

On the MIR/FIR/submm part of Fig. 7 we show the fit to the 
total dust emission, which is the stun of the SEDs from the SF 
regions and the SED from the diflFuse dust obtained in the previ- 
ous sections. The entire dust SED from 8 jum to 850 jum can be 
well fitted. 

On the UV/optical/NIR part of Fig. 7 we show the ob- 
served UV-optical SED (grey squares connected by a solid line). 
Following Eq. C. 12 of Popescu et al. (201 1), the UV/optical SED 
of the young stellar disk of NGC 4214 was derredened using the 
composite attenuation, Am^, which is for the case of old ~ 0: 

Am^ = -2.5 log (1 - F /J + Am/''"* (9) 

where the first part takes into account the attenuation in the SF 
regions and Aw^'"'"* the attenuation of the diffuse component. 
The wavelength dependence of the escape fraction, /^, is tab- 
ulated in Table A.l in Tuff's et al. (2004). For the old stellar 
population we used the attenuation correction derived for the 
old stellar component in the model of Popescu et al. (2011) 
Am^ - Am/'*. For both stellar populations we used the value 
of T-^= 2.0 obtained from our best-fit model and an inclination 
angle of i = 44° (Walter et al. 2008). The SED derived in this 
way represents the total intrinsic stellar SED of NGC 4214 and 
is shown in the blue line. 

We now calculate, in a similar way as in Sect. 6.2 for S FR', 
from the intrinsic UV-to-blue luminosity the predicted value 
of the SFR, by applying Eq. 17 from Popescu et al. (2011), 
SFR ^ L'^'f / lJ^''"."fjy. Using F = 0.40 in Eq. 9, we derive 

SFR = 0.22Moyr \ which is a factor of 2.24 higher than 
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Fig. 6. Fit to the diffuse emission of NGC 4214 keeping fixed old - 0.009 and leaving and SFR' as free parameters. The red solid 
line represents the best fit to the MIR/submm SED for the measured value of the scale-length hs = 873 pc, obtained for t^= 2.0, 
and SFR' = 0.059 M0yr '. The dark grey dashed line corresponds to the sum of the best fit models of the Hll regions presented 
in Figs. 4 and 5 and the grey-filled circles correspond to the sum of the photometric points of NW and SE regions. The light grey 
area represents the uncertainty in the estimation of the total emission of the Hll regions. The value of the reduced (neglecting the 
IRAC 8 fim and MAMBO points) is 1.82. 



the value derived from the fitting of the dust SED (SFR - 
0.098 Moyr '). For F = 0.20 and F = 0.60, the highest and 
lowest values suggested by our data, the corresponding ratios 
are 2.51 and 2.18, respectively. If we adopt for the scale-length 
the maximum value allowed by the data /zj = 1045 pc, we obtain 
the best fit for = 1.2 and SFR' = 0.088 Moyr"', yielding 
SFR = SFR' /(I - F) = 0.088/(1 - 0.40) = 0.14 Mgyr-'. In 
this case, the discrepancy with the value derived from the derre- 
dened UV-optical data, SFR = 0.188 M0yr"', decreases to a 
factor of 1.34 for F = 0.40 (1.50 and 1.27 for the F = 0.20 and 
F - 0.60, respectively). 

A complete discussion of this discrepancy will be given in 
Sect. 7. 

6.4. Gas-to-dust mass ratio 

The total diffuse dust mass can be calculated from Eq. 44 of 



Popescu etal. (2011) (M' 



diff 
dust 



tb X hf X 0.99212 pc^^Mo). With 



our values of /zj - 873 pc and tb = 2 we obtain M^jJ^j = 1 .5 x 
10'' Mq. This value is quite robust again changes in parameters: 
Choosing the largest permitted scale -length, /zs - 1045pc and 



the corresponding best fit for the opacity, Tg - 1 .2, we obtain 
^diilt = 1.3 X lO*" Mo, only 13% lower. From the modelHng 
of the Hll complexes for Z=0.4 Zq, we derived dust masses of 
^dZ = 0.79 X 10^ M© and M^^^ = 0.42 x 10^ M© for the 
NW and SE complexes, respectively. We therefore derived for 
NGC 4214 a total mass of dust Mf"' = Mf«', + M^"^, + Mf = 

dust dust dust dust 

1.62 X 10'' Mq. 

The total gas mass can be derived from Mhi = 4.1 x 10^ Mq 
(Walter et al. 2008) and Mh2 = 5.1 x lO'' Mq (Walter et al. 2001, 
obtained with a Galactic conversion factor) and gives, taking into 
account a helium fraction of 1.36, Mg^j = 5.66 x 10*^ Mq. The 
total gas-to-dust mass ratio is then Gdust = 350. 

The Galactic CO-to-H2 conversion factor (X-factor, adopted 
here as N(H2)//co = 2 x lO^" cm^^ (K km s^i)-'), most 
likely severely underestimates the total molecular gas mass in 
NGC 4214. An indication for this is e.g. the very high Cll/CO 
ratio measured in the central region (Cormier et al. 2010) which 
shows that a large fraction of the CO is photo-dissociated due 
to the high radiation field and low dust shielding. Even though 
the X-factor has been notoriously difficult to determine, progress 
has been made in recent years mainly due to the possibility to 
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Fig. 7. SED for the total emission of NGC 4214. On the MIR/FIR/submm part we show the best-fit solution for the total emission 
(blue solid line) obtained as the sum of the best fit models of the HIl regions (grey dashed line) from Figs. 4 and 5 and the best 
fit model to the diffuse emission (red dashed-dotted line) from Fig. 6, obtained for t^= 2.0 and S FR = 0.098 M0yr '. On the 
UV/optical/NIR part, the grey squares and grey solid line represent the observational data, the blue-filled circles and the blue solid 
line are the intrinsic fluxes derived from deattenuation of the observed data points as described in Sect. 6.3. The orange dotted 
line corresponds to the intrinsic emission of the old stellar population (see Sect. 5.2 for details). The dotted purple line shows the 
intrinsic emission from the young stellar population from the UV to the optical wavelength range. The emission of the young stellar 
population in the B and V band was obtained after subtraction the contribution of the old stellar component (orange dotted line). 
The purple dashed-dotted line is the scaled emission of the young stellar population by the factor needed to predict the same SFR 
as derived from the dust SED modelling. 



derive dust masses with an increasingly better precision. This 
allows, together with HI and CO measurements and the assump- 
tion of a constant gas-to-dust mass ratio, to derive the ratio be- 
tween total molecular gas mass and CO emission. Israel (1997) 
has first used this method and derived for NGC 4214 a 15-30 
times higher X-factor than the Galactic value. Leroy et al. (201 1) 
used spatially resolved data for the dust masses, derived from 
fits to HERSCHEL and SPITZER data, and the gas mass from 
HI and CO measurements for a small sample of nearby galaxies 
(M 31, M 33, LMC, SMC and NGC 6822) . They found evi- 
dence for a strong increase of the X-factor below metallicities 
of 12-log(0/H) = 8.2-8.4, most likely due to the dissociation of 
CO and the creation of extended layers of CO-free H2. Leroy 
et al. (201 1) find for NGC 6822, which has a similar metallicity 
as NGC 4214 (12-Hlog(0/H) = 8.2) an X-factor of 4-5 times the 
Galactic value. For the same galaxy, NGC 6822, and a similar 
method, Gratier et al. (2010) derived a higher X-factor (20 times 
the Galactic value). Taking this range of estimates into account, 
we adopt an X-factor 10 times the Galactic value as a reasonable 
estimate and derive for this case a gas-to-dust mass ratio of 387. 

If we assume that the gas-to-dust mass ratio scales linearly 
with metallicity (which means that the fraction of metals incor- 
porated in the dust is constant) we expect, based on the solar 
value of about 150, a value of Gdust between 375 (for Z = 0.4 Zq) 
to 500 (forZ = 0.3ZQ). 

Thus, the observed gas-to-dust ratio is very close to the ex- 
pected value in the case of Z = 0.4 Z©, and a factor 1.3 lower for 
Z - 0.3 Zq. Even though we consider this discrepancy small, it 
could indicate that there is a still larger amount of "dark" molec- 
ular gas, not probed by the CO. We would need a conversion fac- 



tor 35 times the Galactic value in order to obtain a gas-to-dust 
mass ratio of 500. An alternative possibility is that the dust prop- 
erties in NGC 4214 are different with a higher dust extinction 
coefficient in the submillimeter wavelength range. This could be 
produced e.g. by amorphous graphite as suggested by Galliano 
et al. (201 1) for the case of the LMC. 



7. Discussion 

The analysis presented in this paper enables us to draw conclu- 
sions about the physical properties of the dwarf galaxy under 
study. The way this can be achieved is by confirming or reject- 
ing the basic ingredients of the models used to fit the data, from 
the consistency, between model predictions and data. Since our 
models are radiative transfer models, they contain a wealth of in- 
formation regarding the distribution of stars and dust in galaxies, 
the dumpiness of the ISM, the dust opacity of the dust clouds, 
the contributions of the different phases of the ISM, and the opti- 
cal properties of the dust grains, including the relative abundance 
of PAH molecules. 

Although we can fit the FIR SED of this galaxy, the main dis- 
crepancy is that the UV emission is underpredicted with respect 
to our corresponding predictions for attenuation of UV light. We 
discuss possible causes for this discrepancy, as well as conse- 
quences for the derived physical parameters of this galaxy. We 
also discuss the underprediction of the 8 fim emission in the two 
central star-forming regions, and its implications. 
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7.1. Emission from PAIHs atS/xm 

We found that the observed emission at 8 fim is overpredicted 
for the diffuse emission, and severely underpredicted for the 
emission from the HII regions. One reason for this could be, as 
we discussed in Sect. 3.1, that part of the 8 jjm flux seen in the 
apertures, although powered by UV light from the star cluster, is 
actually likely to originate in the diffuse dust layer beyond the 
PDR+HII region considered by the Groves et al. (2008) model. 
Although it is impossible to quantify the importance of this ef- 
fect, we do not think that it is entirely responsible for the discrep- 
ancies, because most of the 8 jum emission is clearly associated 
with the SF regions, as shown by the close correspondance be- 
tween the 8 /um and Ha distribution (e.g. in the shell structure in 
the NW region.) 

The model of Popescu et al. (2011) assumes Galactic dust 
properties and a fraction of PAHs in dust which is appropriate 
for our Galaxy. The PAH fraction in NGC 4214 might be lower 
given the lower metallicity (e.g. a lower PAH fraction has been 
found in the diffuse ISM in the SMC by Sandstrom et al. 2010). 
If this is the case, we can quantitatively understand the discrep- 
ancy. 

A possible reason for the excess of 8 fim emission in the HII 
regions may be that the model of Groves et al. (2008) overpre- 
dicts the PAH destruction for this particular object. In the models 
of Groves et al. (2008) the PAH abundance is fixed to be pro- 
portional to the gas metallicity. In addition to this, it includes a 
parameter which controls the destruction of the PAHs in intense 
radiation fields. Thus, this model considers PAH destruction and 
assumes a lower fraction of PAHs in dust in the SF regions. 

The observations of Engelbracht et al. (2008) indeed suggest 
that NGC 4214 might be a galaxy with an unusually high 8 yum 
emission for a given radiation field. They found a good correla- 
tion between the equivalent width at 8 jum, EW(8 jum), and the 
ionization parameter for a set of starburst galaxies. The rela- 
tion has also been confirmed for HII regions in MlOl (Gordon 
et al. 2008). NGC 4214, one of the starburst galaxies in the 
Engelbracht et al. (2008) sample, deviates significantly from this 
correlation: it exhibits an EW(8 /im) a factor of 3 higher than 
the one corresponding to the ionization parameter. The spectra 
used in this analysis are not restricted to the SF regions and thus 
might contain emission from the diffuse medium (Engelbracht 
et al. 2008 do not quantify this) and therefore we cannot rigor- 
ously compare this result to our findings for the SF complexes. 
However, since the lines from the ionized atoms and the PAH 
emission come to a large extent from HII regions and their PDRs 
(see our Tabs. 1-3) the result of Engelbracht et al. (2008) sug- 
gests that the 8 fxm emission in the SF regions in NGC 4214 
might indeed be high. 

A final possibility to explain the high 8 fim emission from 
the SF region is the ratio between neutral and ionized PAHs. The 
model of Groves et al. (2008) uses templates fitted to the star- 
burst galaxies NGC 4676 and NGC 7252 with a corresponding 
ratio between neutral and ionized PAHs. This ratio might not be 
entirely appropriate for the low-metalhcity galaxy NGC 4214. 
Ionized PAHs emit about a factor of 10 more energy in the 6- 
9 jum range than at 11-12 fim, whereas neutral PAHs emit about 
the same amount in both ranges (Draine & Li 2007). Thus, a 
higher ionized PAH fraction could boost the predicted emission 
at 8 fim at the expense of the emission at 11-12 fxm. If we as- 
sume as an extreme case that the entire emission at 11-12 fim 
is transferred to the 8 /im band, we can increase the emission 
for our model spectrum at 8 fim by a factor of ~ 2. Thus, we 



could indeed improve the agreement between data and model, 
although some discrepancy would remain. 

7.2. Ttie UV emission illuminating the diffuse dust 

In Section 6 we found that there is tension between the observed 
surface brightness of the disk in the submm and the observed flux 
of the galaxy in the non-ionising UV in the sense that, although 
the FIR SED can be fitted, the UV emission is underpredicted by 
a factor ~ 2-3 with respect to our corresponding predictions for 
attenuation. Here we discuss the plausibility of such an escape 
fraction and study further possibiUties to explain the differences 
between model and data. 

7.2.1 . Escape of unattenuated UV radiation from NGC 4214 

In starburst galaxies, it has been previously observed that a frac- 
tion of the UV luminosity can escape from the galaxies without 
heating the dust, mainly due to the porosity of the ISM created 
by the energy input of massive stars and SNe which in the ex- 
treme case leads to (super)bubbles and galactic winds. Oey et al. 
(2007) showed that starburst galaxies present a low diffuse Ha 
emission consistent with an expected escape fraction of ionizing 
photons of 25% from the galaxy. 

Other studies show that the escape of ionizing radiation (and 
therefore non-ionizing as well) could be much lower: Leitherer 
et al. (1995) pointed out that 3% of the intrinsic LyC photons es- 
cape from a set of starburst galaxies; Grimes et al. (2009) found 
an even lower (<1%) fraction for local starburst galaxies, and 
Siana et al. (2007) found a low escape fraction for low (subso- 
lar) metallicity starburst galaxies (<8%). 

However, the detection of the escaping photons might be de- 
termined by the geometry and orientation of the observations if 
holes and open channels are the medium that the radiation uses 
to escape from the galaxy (Clarke & Oey 2002, Zastrow et al. 
201 1). Besides, the escape fraction can be affected by the inter- 
nal properties of galaxies, depending sensitively on the cover- 
ing factor of clumps, and the density of the clumped and inter- 
clumped medium (Fernandez & Shull 2011). 

Conditions that are favorable for an escape are likely to be 
given in starburst dwarf galaxies which generally show large 
bubbles and shell structures (Martin 1998) - similar to what we 
are seeing in the center of NGC 4214 - that would make the ISM 
more porous and prone to enable the UV radiation to escape. 
The simulations of Wise & Cen (2009) showed that an escape 
fraction as high as 25-80% can be expected due to the irregu- 
lar morphology of the dwarf galaxies with a clumpy ISM. There 
are observational confirmations of such an escape. Cormier et al. 
(2012) found evidence for a escape of 80% of the non-ionizing 
UV-radiation from the low metallicity Haro 11 and Bergvall 
et al. (2006) an escape fraction of the Lyman continuum radi- 
ation of 4-10% in the same object. Zastrow et al. (2011) have 
detected an ionized cone above a shell structure in the starburst 
galaxy NGC 5253 suggesting that ionizing radiation could be 
escaping from this galaxy. In the presence of such channels we 
would expect that not only ionizing, but also non-ionizing radi- 
ation escape from the galaxies without being affected by dust. 

The geometry of the central SF regions, especially in the case 
of the NW region, suggests that bubbles have been already cre- 
ated in NGC 4214. We can make an estimate of how much such 
a process could affect NGC 4214 by assuming that the entire 
UV luminosity from these regions leaves the galaxy unattenu- 
ated. The UV luminosity of the two main HII regions represent 
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20% of the total UV luminosity of NGC 4214. Therefore, about 
20% of the total UV luminosity possibly does not contribute to 
the dust heating and should be subtracted when converting into 
the S FR' heating the dust in the disk. This fraction is close to the 
lower limit of the observed discrepancy and thus might be able 
to explain a considerable part of it. 

7.2.2. A different geometry of dust and stars in NGC 4214 

In our modeling we have assumed that NGC 4214 is a scaled- 
down version of the model spiral galaxy used in Popescu et al. 
(2011). However the differences in the geometry for the spiral 
and dwarf galaxies can affect the comparison of the data with 
the model of Popescu et al. (2011). 

The relative scale-heights and scale-lengths of the dust and 
the stars are derived from a sample of edge-on spiral galaxies 
(Xilouris et al. 1999). However, the geometry of dwarf galax- 
ies is known to be different from that of larger spiral galax- 
ies. In particular, their stellar disks are thicker in comparison 
to their diameter (e.g. Hunter & Elmegreen 2006). This is not 
unexpected because the scale-height depends on the mass sur- 
face density, which determines the gravitational potential in the 
vertical direction and the velocity dispersion. Both are gener- 
ally similar in dwarf galaxies and in spiral galaxies, therefore 
the absolute value of the vertical scale-height is expected to be 
similar, and, since they are smaller than spiral galaxies, the ratio 
between scale-height and scale-length becomes larger. The gas 
scale-height is even higher in absolute terms than in spiral galax- 
ies (Brinks et al. 2002). Thus, we would expect higher ratios be- 
tween the vertical scale-heights and the radial scale-lengths than 
the one assumed in the model. In NGC 4214 Maiz-Apellaniz 
et al. (1999) estimates the vertical stellar scale-height to be about 
200 pc, i.e. about 1/4 to 1/5 of the radial scale-length. This is at 
least a factor of 3 higher than the ratio assumed in the Popescu 
et al. (2011) model. However, the effects of a difference in the 
ratio between scale-height and scale-length are not expected to 
be very large as long as the relative scale ratios between stars 
and dust are the same. 

Another possible, albeit speculative, difference between the 
model and the real geometry could be a non-central position of 
the NW and SE SF regions. Since these regions are dominating 
the young SF activity in the galaxy and are very likely respon- 
sible for a large fraction of the UV emission, their possible po- 
sition above or below the galactic disk would have an important 
impact on the dust heating, whereas in large galaxies such an 
asymmetry is unlikely due to the high number of SF regions. 

The geometry of the galaxy can also affect the way we are 
using the attenuation in the disk. The model for attenuation as- 
sumes that the star-forming complexes are distributed in a thin 
disk, with a radially decreasing exponential distribution. If in- 
stead NGC 4214 had most of the star formation occurring in its 
center, then the model for attenuation in the UV may overpre- 
dict the attenuation if the high local concentration of stars has 
the effect of breaking up the dust layer. 

7.2.3. An extended dust component 

If the dust is more extended compared to the stellar disk than 
what is assumed in the model, such an extended dust component 
will be relatively cold and emit more in the submiUimeter. By 
including such a component we would therefore be able to use 
a higher value of S FR' for the rest of the dust and decrease the 
discrepancy with the observed UV emission. Furthermore, an 



extended, cold dust component could explain the missing flux at 
the long wavelengths. Popescu et al. (2002) have found indica- 
tions of such an extended dust component with ISO data in dwarf 
galaxies, however their data lacked the spatial resolution to di- 
rectly measure its spatial extent. In order to test whether there are 
indications for an extended dust component, we derived the ra- 
dial scale-lengths from SPIRE 250, 350 and 500 jum bands. The 
rather irregular distribution of the dust emission did not allow us 
to define elliptical isophotes at all galactocentric radii, instead 
we derived the median values with the isophotal elliptic annuli 
defined in the g-band image. By carrying out a least square- 
fit to the radial distribution of these median values, we derived 
scale-lengths of 835 pc, 920 pc and 1067 pc for SPIRE 250, 
350 and 500 /im, respectively. The increase of the scale-lengths 
from shorter to longer wavelengths is expected because of the 
decrease of the dust temperature towards the outer parts of the 
galaxy. The longest wavelength (SPIRE 500 //m) best represents 
the scale-length of the dust distribution but most likely still un- 
derestimates it. 

The model of Popescu et al. (20 11) adopts a dust scale-length 
of 1.406 X hg. With hs = 873 pc this would predict a dust scale- 
length of 1227 pc, which is slightly higher than the measured 
scale-length at SPIRE 500 jum band. Thus, although we cannot 
exclude the presence of an extended dust component because the 
dust scale-length can be longer than that at SPIRE 500 jum band, 
we do not find any clear evidence for its presence. 

7.2.4. Different dust properties 

Until now we have mainly discussed different geometrical sce- 
narios for explaining the discrepancies we find in the UV. There 
is however another possibility, that the submm mass absorption 
coefficient of grains, as specified by the Weingartner & Draine 
(2001) model and as used in the radiation transfer model of 
Popescu et al. (2011), underestimates the true mass absorption 
coefficient of the grains in the submm, leading to a lower dust 
surface densities and a lower disk opacity in the UV for the ob- 
served submm brightness of the diffuse disk. This is a solution 
that unfortunately is largely degenerate with the assumed geom- 
etry of stars and dust. Because of this, and the lack of indepen- 
dent constraints on the existence of such a dust model, we cannot 
draw any conclusion on this issue. 

7.2.5. Pliysical implications 

In the previous subsections we discussed a few possibilities that 
might have caused the discrepancy in the predicted level of UV 
emission. Probably the most likely scenario is that dynamical ef- 
fects associated with a recent burst of SF activity responsible for 
the central SF complexes has punched holes through the diffuse 
dust layer, through which the UV photons leaving the SF regions 
can escape the galaxy without interacting with the diffuse dust. 
This would be a natural outcome of a non-steady state behaviour 
of SF activity in dwarf galaxies, in which we may be observ- 
ing NGC 4214 in a post starburst phase, where the star clusters 
are no longer fully cocooned by dust, but to the contrary have 
evolved to the point where SN and wind activity have cleared 
the dust in the ISM above the location of the starburst. 

The existence of different grain properties is also plausible, 
though difficult to distinguish from the other scenarios. This 
would increase the inferred gas-to-dust mass ratios compared 
to the radiation transfer model predictions, bringing it closer to 
the values expected from the observed metalhcity of NGC 4214. 
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Furthermore, dust grains with a higher submillimeter emissivity 
and a flatter frequency dependence would be able to better ex- 
plain the dust emission at long wavelengths (850 /xta) where our 
model presently imderpredict the observations by 34%. 

8. Summary and conclusions 

We have analysed and modelled the dust SED of the nearby star- 
bursting dwarf galaxy NGC 4214. Due to its proximity we were 
able to derive the dust SED separately for the emission from 
the two major massive SF regions and for the diffuse dust. For 
the first time this analysis was done from the perspective of a 
self-consistent radiation transfer calculation constrained by the 
observed SEDs of spatially separated components on resolved 
maps. In making predictions for the UV/optical-FIR/submm 
SEDs of these components this analysis quantitatively takes into 
account the level and colour of UV/optical radiation fields inci- 
dent on both pc-sized dusty structures in the SF complexes and 
on dust grains distributed on kpc scales in the diffuse ISM, il- 
luminated by a combination of UV Ught escaping from the SF 
regions and the ambient optical radiation fields, as constrained 
by UV/optical photometry of the galaxy. The overall analysis 
was done using the model predictions and the formalism from 
Popescu et al. (2011), while for the detailed modelling of the 
two central SF complexes we used the model of Groves et al. 
(2008). The large amount of ancillary data and results from pre- 
vious studies allowed us to constrain a major part of the input 
parameters of the models. 

We achieved a good agreement between data and models, 
both for the diflrise dust emission and the dust in HII-i-PDR re- 
gions. For the comparison with the dust emission from HII-l-PDR 
regions we could constrain from observations basically all in- 
put parameters (metalhcity, age of star cluster, external pressure, 
heating capacity and covering factor) with the exception of the 
gas column density of the PDR. We achieved satisfactory fits 
for both SF regions with the exception of the 8 fim data points. 
Possible reasons for this discrepancy are that the model assump- 
tions (PAH abundance and destruction) are not completely ade- 
quate for the case of NGC 4214 or that NGC 4214 has an unusu- 
ally high emission at 8 /im for its metallicity and radiation field, 
which is supported by other studies (Engelbracht et al. 2008). 

We could fit the diflxise dust SED satisfactorily, but we in- 
ferred that the UV emission was severely underpredicted with 
respect to the observed, deattenuated diffiise UV flux. We have 
discussed diff'erent explanations for this discrepancy (escape of 
UV emission, geometrical effects, a very extended dust disk and 
different dust properties). The most plausible one is that part (30- 
70%) of the UV radiation that escapes from HIl-i-PDR regions 
leaves the galaxy unattenuated and is thus not participating in 
the heating of the diffuse dust, most likely due to a porous ISM. 

We derived a global gas-to-dust mass ratio of 350-390, close 
to the value expected from the metallicity of NGC 4214 (Z = 
0.3- 0.4 Zq) of 375 -500. 

In summary, this is the first time a full radiation transfer anal- 
ysis constrained by the observed SEDs of spatially separated 
components has been done for a dwarf galaxy. Corresponding 
radiation transfer studies of resolved star-forming dwarf galax- 
ies but with different inclinations and evolutionary states to 
NGC 4214 may allow to distinguish between the alternative 
physical scenarios shaping the panchromatic SEDs of these sys- 
tems. 
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